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From molecular spectra to a density shift 
in dense Rydberg gases 
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In Rydberg atoms, at least one electron is excited to a state with a high principal quantum 
number. In an ultracold environment, this low-energy electron can scatter off a ground state 
atom allowing for the formation of a Rydberg molecule consisting of one Rydberg atom and 
several ground state atoms. Here we investigate those Rydberg molecules created by pho- 
toassociation for the spherically symmetric S-states. A step by step increase of the principal 
quantum number up to n = lll enables us to go beyond the previously observed dimer and 
trimer states up to a molecule, where four ground state atoms are bound by one Rydberg 
atom. The increase of bound atoms and the decreasing binding potential per atom with 
principal quantum number results finally in an overlap of spectral lines. The associated 
density-dependent line broadening sets a fundamental limit, for example, for the optical 
thickness per blockade volume in Rydberg quantum optics experiments. 
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Ultr along- range Rydberg molecules, theoretically predicted 
in ref. 1, show very weak binding energies, similar to 
magneto -associated Feshbach molecules^'^, and thus 
require ultracold temperatures. The bond in these molecules 
results from the scattering of a slow Rydberg electron from a 
neutral atom with a negative scattering length a (ref. 1). The 
theoretical approach is based on Fermi's original concept of the 
pseudopotential^ 

Vpseudo (r , R) = <5 (r - R) ( 1 ) 

rriQ 

describing such a scattering event between an electron of mass me 
at position r and an atom at R. If the scattering length a is 
negative, the interaction is attractive and ground state atoms can 
be bound in a potential 

y(R)=^|>P(R)|^ (2) 

me 

where |^(R)| is the local electron density. In case of ^^Rb the 
theoretical value of the scattering length a for triplet scattering is 
— 16.1ao (ref. 5), where ao is the Bohr radius. The singlet 
scattering length is positive and therefore does not lead to a 
bound state. 

The vibrational ground state wavefunction belonging to the 
diatomic (polyatomic) molecules with the highest binding energy 
per atom is localized mainly at the position of the outermost lobe 
of the electron wavefunction, that is, near the classical turning 
point of the electron. Molecules containing more than one 
ground state atom can be described with the same formalism. 
Two ground state atoms inside the electron orbit essentially do 
not interact with each other, as the Rb atom-atom scattering 
length^ is much smaller than the mean particle distance in a 
dilute thermal cloud. Hence, the binding energy of an /-atomic 
molecule, /e N, is (/-I) times larger than the binding energy of a 
dimer. Only when the number of ground state atoms inside the 
Rydberg atom becomes very large, a change in the description of 
the system from discrete bound states to a mean field 
approximation is required. In this case no individual bound 
states are resolved, but the Rydberg line is shifted by 

AE =//drdR Vudo(r,R)|*I'(r)|V(R) 

= /dRy(R)p(R)=2^p, ^> 

where p is the ground state atom density and p its average 
weighted with the probability density of the Rydberg electron. If 
higher order corrections to the zero -energy scattering length 
can be neglected, the mean shift depends only on the value 
of a and p. 

A wide range of fascinating phenomena, starting from the 
coherent creation and breaking of chemical bonds^ to a 
permanent electric dipole moment in a homonuclear molecule^, 
have been already studied in the field of ultr along- range Rydberg 
molecules. Dimers, consisting of a single atom in the Rydberg 
state and one atom in the ground state, have been observed in an 
ultracold gas of Rb in the Rydberg S-state^, D-state^^'^^ and 
P-state^^, and of Cs (ref. 13) in the Rydberg S-state. The regime, 
where the Rydberg electron experiences a mean shift by 
thousands of atoms within its orbit, has been studied at very 
high densities in a Bose-Einstein condensate, leading to electron- 
phonon coupling^ and in a hot vapour^ ^. Here we trace the 
transition between the resolvable molecular lines and the mean 
shift regime in an ultracold cloud with a constant density by 
tuning only one parameter: the principal quantum number n of 
the excited Rydberg state. We probe the limit at which the 
binding energy per atom becomes very small, while the radius of 
the molecules grows beyond 10,000ao. 
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Results 

Molecular spectroscopy. In Fig. 1 excitation spectra from the 
5Si/2 to the nSi/2 state, where n is ranging from 51 to 111, are 
presented. Comparison of some of the relevant molecular 
potentials are shown in Fig. 2. The shape of the obtained Rydberg 
spectra varies significantly for different n. For low principal 
quantum numbers clearly distinguishable molecular lines are 
present on the red side of the atomic peak, which is situated at the 
origin. We can identif)^ the observed molecular lines by 
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Figure 1 | Overview of the 5Si/2 to nSy2 excitation spectra. The dimer— 
polyatomic molecule transition for increasing principal quantum number is 
visible. The origin of the relative frequency axis corresponds to the centre of 
the atomic line (dashed line) for n<7^, where molecular lines are 
distinguishable. The grey shaded area between —30 and 30 kHz indicates 
the laser bandwidth. Spectra at n>7^ are horizontally shifted such that their 
centres of gravity overlap with the mean centre of gravity <cg> (orange 
dashed line) of the first three spectra (n = 51, 62, 71). All data were taken at 
similar cloud parameters, therefore the density-induced shift for all spectra 
is constant in the first approximation. Molecules with up to three bound 
ground state atoms for 625 and up to four for 715 are resolvable in the 
spectra. Coloured diamonds indicate the positions of the dimers (red), 
trimers (violet) etc. following the power law scaling of the binding energies 
fitted to the first three spectra. In the inset, the molecular spectrum for the 
n = 71 Rydberg state is shown. A multilorentzian fit (green line), assuming a 
constant spacing between the molecular peaks is plotted to indicate the 
positions of the higher order molecular lines. The spectrum for n = 40 is not 
shown, because the binding energy of the dimer is larger than the 
plotting range. Each spectrum is an average over 20 independent 
measurements with standard deviation error bars. arb. u., arbitrary unit. 
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Figure 2 | Scattering potentials for n = 51, 62, 71 and 80 5-states. The 

potentials are calculated from equation (2). Dimers are predominantly 
bound in the vibrational ground state in the outermost well (probability 
density in light blue). At n = 80, the mean binding energy per atom (light 
blue dashed line) becomes comparable to the experimental resolution and a 
single-molecular state localized in the outermost lobe cannot be resolved in 
the experiment. The spheres below illustrate the highest order polyatomic 
molecule observed for each principal quantum number. The orange dashed 
line indicates the mean shift for the experimental density. All potentials, 
binding energies and probability densities are plotted to scale. 

calculating the binding energies of dimers in the ground and the 
first vibrationally excited state^. The binding energy can be 
directly measured as the difference between the atomic and the 
molecular line in the spectrum. In the spectrum of n — 5\ the 
peak at — 1.7 MHz can be identified as a dimer, for which the 
ground state atom is localized in the outermost well of the 
molecular potential. Additionally, at —1.3 MHz an excited 
vibrational state^'^^ is visible. For a given density, the 
probability to find an atom inside the electron orbit scales as 
n*^ with the effective principal quantum number rf — n — Sq. 
Hence, higher order molecules are formed more likely at higher n. 
At the same time the binding energy per atom decreases and thus 
polyatomic molecular lines become visible in the spectra of n = 62 
and 71. At n — 62, lines up to the tetramer together with 
corresponding vibrational excited states are visible. The 
broadening of the tetramer line may be caused by the presence 
of these excited states, possibly with a reduced lifetime^^. At 
n = 71 only vibrational ground states are resolved. Polyatomic 
molecules up to a pentamer can be identified. The size of such a 
molecule becomes enormous due to the Rydberg electron orbit 
radius reaching almost 10,000ao (Fig. 1). Since the thermal cloud 
is much bigger than the Rydberg atom, the Gaussian density 
distribution does not play a role on a single Rydberg atom level, 
and the spatial arrangement of the bound ground state atoms is 
isotropic. For large n the binding energy decreases until it is 
below the experimental resolution. This manifests in a non- 
resolvable shoulder and finally in an inhomogeneously broadened 
spectral line. 

Analysis of the spectra. The experimental binding energies can 
be calculated based on the molecular potential (2) and the mean 
shift (3). Momentum-dependent corrections to these values can 
be estimated using a semiclassical approximation^^. Additionally, 
a p-wave shape resonance can cause a substantial modification of 
the molecular potential, leading to butterfly-type molecules^^ and 
states bound by internal quantum reflection . These corrections 
are largely negligible for high principal quantum numbers and 
large distances from the Rydberg core, where the relative motion 



of the Rydberg electron and the perturber is slow. We use 
corrections to the s-wave scattering length including terms linear 
in the relative momentum. As this approximation is valid for 
large distances from the ionic core, we restrict the analysis to the 
lowest bound state. The discussion of higher vibrational states can 
be found in ref. 16. We solve the Schrodinger equation for the 
ground state atom in the molecular potential using Numerov's 
method to find the binding energies of the dimers. We fit the 
zero-energy scattering length a to obtain the best agreement with 
the measured values. In this paper we achieve the best agreement 
with the experimental data for — 16.2ao, which is very close to the 
theoretically predicted value of — 16.1ao (ref. 5). 

For n>71, where no distinct molecular lines can be identified, 
a mean field description is required. Furthermore, the spectral 
position of pure Rydberg atoms, and thus the zero position, 
cannot be identified directly from the signal. However it can be 
determined from the centre of gravity eg of the spectra, taken in a 
non-blockaded sample, which is constant for a given density. 
Intuitively, this result in the first approximation can be explained 

by the fact that while the mean potential depth V — /^^ ^^^-^ 

J 

averaged over the volume of the Rydberg atom decreases with the 
effective principal number as the probability to find an 

atom inside the Rydberg electron orbit increases with n*^. In the 
experiment all data were taken at a fixed density. Thus, in our 
analysis we overlap the centre of gravity of the spectra at n>71 
with the mean <c^> = —300 kHz determined at low principal 
quantum numbers. Doing so we can identif)^ the zero position in 
the top panels of Fig. 1. Assuming the scattering length a to be 
constant, we determine the effective density to be 3 x 10^^ cm~ ^, 
which is close to the peak density obtained from a Gaussian fit to 
absorption images of the thermal cloud. This indicates that 
molecules are most likely created in regions of high density. Note 
that the power broadening would not effect the centre of gravity 
of the measured spectra. Only for n<80 the highest signal 
originates from pure Rydberg atoms. Already for n = 80 there is 
on average one ground state atom inside the electron orbit, 
leading to a high probability to excite dimers instead of pure 
Rydberg atoms. On average there are four atoms inside the lOOS 
electron orbit and eight for the 11 IS state. Therefore, the atomic 
line is suppressed, while the molecular contribution to the overall 
signal increases. At the same time the molecular lines are not 
resolvable in the experiment any more, because of their very low 
binding energies. For even higher principal quantum numbers the 
shape of the spectrum is expected to become Gaussian, with the 
maximum at the position of the mean shift (<c^) in Fig. 1) 
(refs 14,20). 

The binding energies of all observed molecules are plotted in 
Fig. 3. From the extrapolated molecular binding energies the 
transition from the few particle description of discrete bound 
states to a mean field shift becomes visible. A power law fitted to 
the n>40 data shows a scaling with the effective principal 
quantum number n* to the power of — 6.26 ±0.12, close to the 
value of — 6, expected from the size scaling argument. The 
deviation can be explained by the dependence of the scattering 
length on the relative momentum and the fact that with 
increasing principal quantum number the shape of the outermost 
well of the molecular potential changes. Taking the corresponding 
zero-point energy and momentum- dependent corrections to the 
scattering length into account, we obtain an exponent of — 6.37, 
which is in very good agreement with our experimental data. 
Contributions of higher order partial waves and p-wave shape 
resonances to the molecular potential can be neglected, as the 
kinetic energy associated with the relative motion of two 
scattering partners is small in the region of interest. 
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Figure 3 | Measured binding energies versus principal quantum number 
and diameters of the molecules. The data for n<37 are taken from ref. 9. 
For n = 40, 43, 51, the frequency range chosen in the experiment was too 
small to photoassociate molecules with larger binding energies than a 
dimer. The power law n*^ (blue lines) is fitted to the measured binding 
energies (red points) of the dimers and multiplied by factors /-I, ie N, (/ = 2 
for a dimer). The grey shaded area indicates the mean shift regime, which 
starts at around n = 75, which can also be seen in Fig. 1. At this principal 
quantum number, the binding energy of the dimer becomes smaller 
than the experimental resolution. Calculated values of Eq are not shown in 
the plot, because they are hardly distinguishable from the experimental data 
on this scale. The orange dashed line indicates the calculated {eg} from 
the measured spectra (Fig. 1). The error bars are determined as the 
standard deviation of the fit. 



Discussion 

The obtained experimental results have strong implications, for 
example, in the field of Rydberg quantum optics. The density- 
dependent dephasing resulting from the existence of many 
molecular lines within the Rydberg line envelope sets a 
fundamental limit for the number of atoms inside the blockade 
radius, that is, the optical thickness, for n>80. This fact is of 
importance for every experiment taken at high principal quantum 
numbers and high densities, in particular for quantum optics 
experiments in ultracold clouds^ ^ . In many experiments based 
on electromagnetically induced transparency using a Rydberg 
state, the figure of merit is high optical depth per Rydberg 
blockade volume, which is set by the width of the 
electromagnetically induced transparency window^^. Assuming 
a typical linewidth of 1 MHz, density- induced decoherence 
becomes relevant at the density lO^^cm"^ (ref. 23). Therefore, 
this effect has to be taken into account in the experimental work 
and proposals, as in practice it will limit the Rydberg blockade. 

Nonetheless, Rydberg molecules offer a unique tool to address 
few-body subsystems with control on a single-particle level by 
changing the detuning of the excitation laser light. The number of 
constituents of a molecule and the interaction strength can also be 
varied independently by changing the density and the principal 
quantum number of the excited Rydberg atom. Furthermore, the 
analysis of the relative strength of the molecular lines opens up 
the possibility to measure correlations in a bosonic gas. In high- 
density gases and for low principal quantum numbers, where the 
size of the Rydberg atom is comparable to the de Broglie 
wavelength, extracting the ^^^^ correlation function^^ of thermal 
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and Bose-condensed gases is feasible. In addition to previous 
measurements^ ^'^^ also higher order correlation functions can be 
studied using polyatomic molecules. 

Metliods 

Experiment. We perform the experiment in a magnetically trapped ultracold cloud 
of ^^Rb atoms in the 5Si/2, F = 2, mp = 2 ground state with typical temperatures of 
2 |j.K and densities on the order of lO^^cm"^. Detailed information about the 
setup can be found in ref. 29. After the preparation of the ultracold cloud, we excite 
the atoms in a two-photon process via the 6P3/2 state to the nSi/2 Rydberg state 
(quantum defect ^0 = 3.1311807(8) (ref 30)). The detuning from the intermediate 
state in the two -photon transition is 80 MHz. The light driving the lower transition 
is sent to the experiment in 50 |is pulses with a repetition rate of 167 Hz. During the 
sequence the light driving the upper transition is on constantly. Both laser beams 
have beam waists of 0.5 mm. After the excitation we field-ionize the Rydberg atoms 
and collect the ions on a microchannel plate detector. In a single atomic cloud we 
perform typically 400 cycles of Rydberg excitation and detection while scanning the 
frequency of the infrared laser light. We choose long excitation pulses of 50 |is to 
not be limited by Fourier broadening. Taking further into account the Doppler 
broadening this results in an experimental resolution with full-width at half- 
maximum value of 100 kHz. To obtain the best visibility while changing the 
principal quantum number of the excited Rydberg state, we adjust the power of the 
laser driving the lower transition to account for power broadening, therefore we 
maintain a two-photon Rabi frequency Q below 600 Hz. Only the spectrum of 
n = 51 was taken with higher laser power (Q » 750 Hz) and thus in this case the 
atomic line is slightly broadened. With these Rabi frequencies the calculated 
blockade radii for the measured states are ranging from 7|im (5 IS) to 28|am 
(11 IS); for comparison the 3a radius of the thermal cloud is around 0.4 mm. All 
the measurements were taken in the dilute regime for Rydberg atoms. 
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